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Abstract 

Envelope glycoproteins of Hepatitis C Virus (HCV) play an important role in the 
virus assembly and initial entry into host cells. Conserved charged residues of 
the E2 transmembrane (TM) domain were shown to be responsible for the 
heterodimerization with envelope glycoprotein E1. Despite intensive research 
on both envelope glycoproteins, the structural information is still not fully 
understood. Recent findings have revealed that the stem (ST) region of E2 also 
functions in the initial stage of the viral life cycle. We have previously shown the 
effect of the conserved charged residues on the TM helix monomer of E2. 
Here, we extended the model of the TM domain by adding the adjacent ST 
segment. Explicit molecular dynamics simulations were performed for the E2 
amphiphilic segment of the ST region connected to the putative TM domain 
(residues 683-746). Structural conformation and behavior are studied and 
compared with the nuclear magnetic resonance (NMR)-derived segment of E2 
(2KQZ.pdb). We observed that the central helix of the ST region (residues 689 
- 703) remained stable as a helix in-plane to the lipid bilayer. Furthermore, the 
TM domain appeared to provide minimal contribution to the structural stability 
of the amphipathic region. This study also provides insight into the orientation 
and positional preferences of the ST segment with respect to the membrane 
lipid-water interface. 
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BHgSSHD Changes from Version 1 

The following revisions have been made to the article in response 
to the reviewer comments: two new sections added to the 
manuscript "Residue contributions and relevance" and "Limitations 
and future work", the term "membrane interface" has been 
replaced with "lipid-water interface", the conclusion has been 
rephrased and the first structure of the NMR ensemble has been 
added to the "Input structure preparation" section. 

See referee reports 



Introduction 

Envelope glycoproteins El and E2 are essential for the initial bind- 
ing and internalization of Hepatitis C Virus (HCV) into host cells. 
Both glycoproteins have been shown to interact as a non-covalent 
heterodimer during biosynthesis 1 . Several conserved charged resi- 
dues located in the transmembrane (TM) domains of El and E2 
were shown to function not only as membrane anchors, but were 
also essential for dimerization, endoplasmic reticulum retention and 
viral envelope formation 2 ' 3 . 

In our previous studies, we demonstrated that the unfolding be- 
havior of the E2 TM helix monomer was attributed to the charged 
Asp728, which was located in the hydrophobic core. The main 
contribution of Asp was postulated to be located at the helix-helix 
interface and involved the formation of a salt bridge with the Lys 
of the El envelope glycoprotein 2 . The ion-pair interaction of the 
E1-E2 heterodimer was captured in the molecular dynamic (MD) 
simulation studies based on the model that placed the charged Asp 
and Lys at the helix-helix interface 4 5 . 

E2 envelope glycoprotein is known to be required for interactions with 
cellular receptors involved in endocytosis and membrane fusion 6-8 . 
E2 is composed of domain I-III, followed by the ST region and 
the TM domain 7 . Recent studies based on circular dichroism (CD) 
and nuclear magnetic resonance (NMR) spectroscopy revealed that 
the soluble region located adjacent to the TM domain of E2 was 
involved in the initial virus entry. This highly conserved ST region 
was shown to fold as a helix upon membrane binding 9 . 

In this work, we carried out MD simulations for three E2 structures: 
(1) a model generated by the I-Tasser server, (2) an ideal helix model 
and (3) an NMR derived structure of the E2 segment (2KZQ.pdb). 
The first two models include the TM domain of E2 but the TM do- 
main is not present in the NMR structure. We observed consistent 
structural stability in the ST region amphipathic segment (residue 
689-703) across all simulations suggesting that the contribution of 
the TM domain to the segment structure stability is minimal. In ad- 
dition, we demonstrated the orientation and positional preferences 
of this amphipathic segment. 

Methods 

Input structure preparation 

The protein sequence of HCV E2 genotype la (H77 strain) (Uniprot 
ID P27958) used to prepare the models for MD simulations was ob- 
tained from UniProtKB/Swiss-Prot database (www.uniprot.org) 10 . 



The following is the sequence segment of HCV E2 used to prepare 
the first model for this work, referred to as an ideal helix model: 
683 PALSTGL 690 IHLHQNIVDV 700 QYLYGVGSSI 710 ASWAIKW- 
EYV 720 VLLFLLLADA 730 RVCSCLWMML 740 LISQAEA. The ideal 
helix model was generated using Pymol (http://pymol.sourceforge.net) 
by orienting the ST segment (residue 683-714) perpendicular to the 
TM domain (residues 715-746). The protein structure images in this 
work were prepared by the Pymol program. The I-Tasser webserver 11 
was used to generate the second model. The complete sequence of 
E2 was submitted to the I-Tasser server. Five models were generated 
and the model with a low root mean square deviation (RMSD) with 
the available NMR structure (2KZQ.pdb) and consisting of a helical 
TM domain was selected. Only the same segment as examined with 
the ideal-helix model was used for further MD simulations. The first 
model of the E2 NMR ensemble (PDBID: 2KZQ) retrieved from the 
Protein Databank (PDB) 12 was used in the third simulation system. 
This E2 protein segment was based on the HCV genotype 2a (JFH-1) 
(Uniprot ID Q99IB8) 9 . 

System preparation 

Pre-equilibrated dipalmitoylphosphatidylcholine (DPPC) lipid bilayer 
was retrieved from the web of Prof. Tieleman (http://moose.bio. 
ucalgary.ca/). Peptide orientation in DPPC lipid bilayer was done by 
aligning the hydrophobic belt of the peptide, parallel to the membrane 
plane using LAMBADA 13 . The optimal number of overlapping lipid 
molecules was subsequently calculated and removed followed by 
lipid expansion (inflation) and alternating twenty steps of deflation and 
energy minimization to allow the peptide to be embedded within the 
bilayer using inflateGR02 13 . A short 100 ps energy minimization was 
employed to relax possible steric conflicts. Ions and counter ions were 
added to neutralize the system followed by 20 ns position-restrained 
simulation allowing the bilayer to re-equilibrate around the protein. 
Production MD simulations were carried out for 100 ns in the I-Tasser 
model, and 20 ns for both the ideal helix model and NMR structure 
(2KZQ.pdb). 

MD simulations 

The DPPC lipid bilayer interactions were described using the 
Berger force-field parameters 14 . The TM helices were modeled with 
the united atom force-field GROMOS96 53a6 15 . Simulations were 
performed with the Gromacs 4.5.5 package 16 using 2-fs time steps. 
Periodic boundary conditions were used in all directions. Bonds to 
hydrogen atoms were constrained using the LINCS algorithms 17 . 
For the short-range van der Waals interactions, a cut-off distance 
of 1.0 nm was used. The long-range electrostatic interactions were 
treated using the particle mesh Ewald (PME) method with a grid 
spacing of 0.12 nm and cubic interpolation. The non-bonded pair 
list was generated every 10 steps with a cut-off of 1.0 nm. Water, 
lipid and peptide systems were coupled separately to temperature 
baths with 323 K for the DPPC using the Berendsen algorithm with 
a time constant of T x = 0.1 ps 18 . To maintain constant pressure, semi- 
isotropic coupling was employed separately for the lateral and for 
the normal directions with Berendsen weak coupling and a T p = 1 
ps time constant. The compressibility was set to 4.5 x 10~ 5 bar 118 . 

Analyses of the trajectories were primarily performed with tools 
included in the Gromacs 4.5.5 suite 16 . RMSD analyses were based 
on the coordinates of all atoms of the peptides. The bilayer thickness 
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was measured by averaging the distances between lipid headgroups 
in the upper and lower leaflets of the lipid membrane with the tool 
GridMAT-MD 9 . 

Results and discussion 

Stability of the amphipathic segment 

A stable helical conformation of the E2 ST region (residue 689-703) 
was consistently observed with some uncharacteristic spikes in the 
early stages and towards the end of the MD simulation of the I-Tasser 
model and the NMR derived structure (2KZQ.pdb), respectively. 
RMSD of this amphipathic segment was also consistently observed to 
be progressing within the commonly accepted 2 A range for the ideal 
helix and NMR structures throughout the simulations. On the other 
hand, the I-Tasser model showed subtly higher RMSD progression 
over the simulation time (Figure 1, Figure 2 and Figure 3 a, 3b). 

Secondary structure stability observed in these three contrasting 
simulation systems can be attributed to the amphipathic nature 
of the residues allowing the helix to retain its structure on both 
the hydrophobic core of the lipid bilayer and the hydrophilic 
environment of the solvent. Furthermore, secondary structure 
analyses done by Albecka et al. predicted that the amino acid 
segment 689-703 consists of more than 50% helix. In addition 
to that, consistent results were shown from CD spectra and NMR 
with 50% 2,2,2-trifluoroethanol (TFE) environment, which both 
detected the presence of helical segments in the E2-ST region. 
The most amphipathic helix characteristics are displayed in the 
E2-ST segment 687-703. Based on the profile, Albecka et al. 
speculated that this segment is a helix upon membrane binding 9 . 
Interestingly, the presence of the TM domain in our simulation 
does not appear to contribute significantly to the helix stability of 
the amphipathic region. 

The helical length of this region does not vary much between the 
I-Tasser and the ideal helix models (Figure 2 and Figure 3a), which 
include both the ST region and TM domain, compared with the 
2KZQ.pdb (Figure 3b) which does not include the TM domain. 
This observation suggests that lipid-peptide interactions play a 
larger role in stabilizing the secondary structure of this amphip- 
athic segment compared with the TM domain. These data provide 
evidence of the contribution of lipid to structural stability modula- 
tion and are in good agreement with the hypothesized lipid and/or 
protein contribution to structural stability 9 . The higher RMSD value 
observed in the I-Tasser model simulation was well anticipated and 
was mainly attributed to the relative positioning of the ST region, 
which was sandwiched in between lipid leaflets, forcing the seg- 
ment to reorganize its structure conformation and having only a 
minimal effect on the helical integrity of the secondary structure. 
Examining this reorganization further by monitoring the distance of 
the amphipathic segment to lipid leaflets led to another interesting 
observation described in the next section of this article. 

Orientation and positional preference of the amphipathic 
segment 

Monitoring the movement of the amphipathic segment during 
the simulations led to another interesting observation. Segment 
of residues 689-703 in the I-Tasser model appeared to move to- 
wards the hydrophobic core of the lipid bilayer as depicted by 



Figure 1. Superimposition of the amphipathic region of E2 
(segment 689 to 703) in three molecular dynamics simulations. 

The most stable helical region of the E2 amphipathic segment 
during molecular dynamics simulations. 2KZQ in red, I-Tasser model 
in black and ideal helix model in gray. 
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Figure 2. Helix length and RMSD of I-Tasser model. Both 
amphipathic regions and transmembrane domain located in the 
hydrophobic core of the bilayer. Length of helix (residues 689-703) 
is plotted in red, root mean square deviation of the same residues 
with respect to the starting structure is plotted in black. 

steady progression in the distance to both the upper and lower 
lipid leaflets depicted in Figure 4. This reorganization is surpris- 
ing and interesting because we would have previously assumed 
that the amphipathic region exposed to the solvent would hold 
the structure steady, despite some part of the amphipathic region 
being initially vertically positioned in the hydrophobic core of the 
lipid (Supplementary Figure 1). In addition, given the amphipathic 
nature of the residues in this segment, one could postulate that the 
residues would remain in this position. However, over the period 
of the simulation the segment reoriented by moving away from the 
lipid leaflets, while at the same time retaining structural integrity. 
The structural stability of the segment, as discussed in the previ- 
ous section of this article, is attributed to the amphipathic nature 
of the residues but this does not explain the movement towards 
the hydrophobic core of the lipid. The systematically orchestrated 
movement towards the hydrophobic core of the lipid leaflets indi- 
cates a strong orientation preference of the amphipathic segment, 
which in this specific case was parallel to the lipid leaflets. We 
then monitored the segment movement relative to the lipid leaf- 
lets with the other two simulations (the ideal helix model and the 
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Figure 3. Helix length and RMSD of the ideal helix model and 2KZQ structures, (a) Molecular dynamics simulation of the ideal helix 
model. The amphipathic segment was positioned in water-lipid interface and the transmembrane domain was oriented in the hydrophobic 
core of the bilayer. The length of the helix (segment residue 689-703) is plotted in red, root mean square deviation (RMSD) of the same 
residues with respect to the starting structure is plotted in black, (b) Molecular dynamic simulation of the nuclear magnetic resonance derived 
structure (2KZQ.pdb). 2KZQ was positioned in the water-lipid interface. The length of the helix segment is plotted in red; RMSD of the same 
residue with respect to starting structure is plotted in black. 



2KZQ structure). The segment was initially positioned horizontally 
to the lipid leaflets. The results showed that the distance of the 
amphipathic segment in both simulations was consistently within 
4 A to the lipid phosphate head group throughout the simulation 
(Figure 5). These data further clarify the orientation preference 
of the amphipathic segment with respect to the lipid leaflets and 
suggest that the residues are positioned at the lipid- water interface 
in a very stable manner. Interestingly, Albecka et al. 9 speculated 
that these residues could have an in-plane topology or orientation 
and suggested that the ST region would ideally be positioned at 
the membrane lipid-water interface, which is again in agreement 
with our data 9 . 

Residue contributions and relevance 

Having observed a clear preference of amphipathic segment ori- 
entation in the lipid bilayer environment, it was also pointed out 
to us that it would be desirable and interesting to look into residue 
contribution on this behavior. We then plotted the residue distribu- 
tion in the lipid water interface for the 2KZQ structure and ideal 
helix model (Supplementary Figure 2) to identify specific residues 
positional preference in the lipid water interface. As depicted in 
Supplementary Figure 2, a clear residue preference in lipid-water 
interface was not observed. We found polar and charged residues 
distributed in both hydrophilic and hydrophobic regions of the in- 
terface; a similar observation was also found for non-polar residues. 
This result indicates the overall contribution of the residues played 
more roles in the stability and positional preference of this segment 
in lipid-water interface superseding any individual residue contribu- 
tion. Characterizing the residue distribution in our I-Tasser model 
would not provide information on the residues preferences since the 
segment is situated exclusively in a hydrophobic environment. To 
get some insight on this segment, residue rotation analysis of this 
segment coupled with qualitative visualization of polar, nonpolar 
and charged residues was carried out (Supplementary Figure 3). 
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Figure 4. Distance between the helix (residues 689-703) and 
lipid leaflets during molecular dynamic simulation of the 
l-Tasser model. The distance to the upper lipid leaflet is plotted in 
red and lower leaflet is plotted in black. 

Synchronized rotation of the amphipathic segment was observed 
and convergence was captured after 50 ns and for the remaining 
period of the simulation. In contrast to the ideal helix model and 
2KZQ structure simulations, a clear preference to hydrophilic and 
hydrophobic regions was observed. Polar and charged residues 
were oriented towards the hydrophilic region and non-polar residues 
were oriented towards the hydrophobic region. We also noticed that 
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Figure 5. Distance of the amphipathic segment to the lipid 
leaflets. The 2KZQ is colored in black while the ideal helix model 
is colored in red. Amphipathic residues remain within 4 A from the 
phosphate head group throughout the simulations. 



convergence in rotation on 50 ns (Supplementary Figure 3) coincide 
nicely with a major shift in distance profile (Figure 4), which also 
started at around 50 ns and stabilized at around 70 ns. This indi- 
cates that the reorganization of this segment correlates tightly to the 
residues' preference towards hydrophilic or hydrophobic environ- 
ments. It is interesting to highlight that these residues' preference to 
hydrophilic and hydrophobic regions in lipid-water interface (ideal 
helix model and 2KZQ structure) appears to be loosely correlated 
(no clear preference observed) in contrast to the residues' prefer- 
ence in environments exclusively on lipid (I-Tasser model) where 
a clear preference was observed. Due to the fact that the movement 
of the I-Tasser model residues is synchronized, we found no evi- 
dence to substantiate major contribution of individual residues to 
this structural rearrangement. Overall hydrophilic and hydrophobic 
residue contribution to the segment reorientation and reorganiza- 
tion is further highlighted. 

Limitations and future work 

An orientational preference parallel to the lipid leaflets is a com- 
mon feature of amphipathic helix segments. Parallel orientation of 
amphipathic segments has been characterized in experimental 
works 20 . Amphipathic segments have also been observed to only 
minimally penetrate and contribute to lipid bilayer perturbation 2122 . 
In addition to that, lipid shape has been demonstrated to contribute 
to amphipathic segment topology in the lipid interface 23 . Taking all 
this information collectively, the amphipathic segment of residues 
689-703 is anticipated to be positioned stably in the lipid-water in- 
terface in a parallel topology with respect to the lipid leaflets. In line 
with this anticipation, the 20 ns simulations of two of our structures 
(ideal helix model and 2KZQ) demonstrate a consistent helical form 
oriented parallel to the lipid leaflet throughout the simulation period. 



Interestingly, our third simulated structure (I-Tasser model) demon- 
strated a clear preference to parallel orientation despite the fact that 
its initial amphipathic segment configuration was perpendicular to 
the lipid leaflet. We decided to subject the I-Tasser model to longer 
simulation (100 ns) compared to the other two structures (20 ns) 
mainly due to the position of its amphipathic segment. The I-Tasser 
model positioned the amphipathic segment sandwiched in the bi- 
layer leaflets, in contrast to the commonly assumed position in the 
membrane interface. On the other hand in the ideal helix and 2KZQ 
models it is readily positioned at the membrane interface. True to our 
expectation, meaningful observations such as amphipathic segment 
reorganization only started to be captured after 70 ns. This highlights 
the importance of optimum conformational space exploration in the 
simulation. Despite the consistent stability observed throughout our 
ideal helix model and 2KZQ structures, a longer simulation time is 
desirable to further reaffirm the observation. Simulation of the am- 
phipathic segment in the water-only environment is also desirable to 
elucidate the segment behavior. These will be added as part of the 
future agenda of this work. 

Conclusion 

In this study, the atomistic MD simulations provide insightful struc- 
tural data for the E2 segment. The amphipathic segment of E2 was 
able to remain as a stable helix in a lipid bilayer environment even 
without the respective TM domain. The results also describe ori- 
entation and positional preferences of the amphipathic segment in 
relation to the lipid-water interface. The amphipathic segment con- 
sisting of residues 689-703 demonstrated distinct parallel positional 
preference to lipid-water interface. Dominant individual residue 
contribution was not observed in the segment rearrangement and 
reorganization, overall residue preference in this rearrangement and 
reorganization is highlighted. 
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Supplementary figures 




Supplementary Figure 1. Molecular dynamic simulation of l-Tasser model, (a) Initial starting (cyan) and (b) 100 ns snapshots (green) of 
molecular dynamic simulations. Segment residues 689-703 are colored in red. 
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Supplementary Figure 3. Residue rotation of l-Tasser model, (a) Synchronized rotation was observed, convergence was also observed at 
around 50 ns and the remaining of the simulation time, (b) Residues pointing to hydrophilic regions are denoted with asterisks and colored in 
cyan and residues pointing to hydrophobic regions are colored in yellow. 
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Approved: 03 September 2013 
Referee Report: 03 September 201 3 

I appreciate the effort by the authors to improve their draft following the reviewers suggestions. Although 
they did not address all the points mentioned, I think the current version is more consistent, and few 
vague comments have been modified. The additional sections provide useful information on the system 
that improve the understanding of the preference in position and orientation of the amphipathic segment. 

Regarding the extension of the simulations, although I would keep insisting in the importance of it, I can 
understand that the required computer power might not be accessible to everyone. They have included a 
sentence pointing at this factor, and as far as I am concerned, it should be enough for the acceptance of 
their work. I am looking forward to the extension in this topic from the authors. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 

Competing Interests: No competing interests were disclosed. 




Martin Zacharias 

Physics Department, Technical University Munich, Garching bei Munchen, Germany 



Approved: 12 August 2013 
Referee Report: 12 August 2013 

I believe the authors made additions and corrections to the original article following the suggestions of the 
first review report. The article is now acceptable including a discussion on the limitations of MD 
simulations length and force field approximations. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 

Competing Interests: No competing interests were disclosed. 
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2 CNB-CSIC, Madrid, Spain 

Approved with reservations: 18 April 2013 
Referee Report: 18 April 2013 

General Comments and Recommendation: The title is appropriate, the design of the simulations, and 
the methods and analysis employed are state of the art. The data provided is sufficient. However, the 
conclusions are overstretched and much longer simulations are needed to substantiate the claims. 

More detailed review:The article describes simulation studies of a 64-residue transmembrane fragment 
of the envelope glycoprotein E2 domain consisting of two segments: TM (residues 71 5-746) and ST 
(residues 683-714). The objective of the study is mainly to understand the role and structural dynamics of 
the ST segment. MD simulations have been performed on three different starting points: 2 models and 
one experimentally derived structure. The experimental structure 2KZQ is only 36 residues long itself, 
lacking the TM segment, and has large dynamics within the NMR structural ensemble (36 models with 
mean RMSD of 1 1 .32). It is the homolog that produced the best model in l-Tasser, the authors have used 
in the study. The authors do not mention which model from the NMR ensemble they used in their MD 
simulations. Though the authors refer to previous experimental studies that show that the ST segment is 
helical upon membrane binding, a simple sequence-based helical propensity prediction of the HCV E2 
genotype 1a (H77 strain) sequence (Uniprot ID P27958) could be included (to further substantiate their 
'ideal helix 1 model used for the MD simulations). One of the main results of the article is that the ST 
segment has a stable helical conformation, stabilized more by the lipid-peptide interactions than the 
TM-ST peptide-peptide interactions. The similarity between the three different simulations, even with one 
of them lacking the TM segment, forms the basis of their conclusion. The authors observe 
uncharacteristic spikes towards the end in the simulations of NMR-structure and the l-Tasser model MD 
simulations (Figures 2, 3a, 3b). Since the NMR-structure simulation is very short (just 20ns) compared to 
the l-Tasser model simulation, the authors could perform a longer, equivalent 100ns simulation to 
compare them and ascertain further the importance of lipid-peptide interactions in the ST helical segment 
stability. Also, given that in the l-Tasser simulation, there is a lot of dynamics observed after the 20ns 
timepoint until the end of the simulation, a longer simulation of the NMR-structure is recommended. The 
other key result which we found interesting is the orientational and positional preference of the 
amphipathic ST segment, which seems to be parallel to the lipid leaflets. Reorganization and movement 
of the ST segment towards the hydrophobic core of the lipids, as shown in Fig S1 , is observed clearly in 
the l-Tasser model simulation, whereas the initial horizontal orientation is maintained in the other 
simulations. The authors could identify the key residues involved in this re-organization and highlight them 

We have read this submission. We believe that we have an appropriate level of expertise to 
confirm that it is of an acceptable scientific standard, however we have significant reservations, 
as outlined above. 

Competing Interests: No competing interests were disclosed. 
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Gloria Fuentes 

Bioinformatics Institute, A*STAR, Matrix, Singapore 

Approved with reservations: 18 April 2013 
Referee Report: 18 April 2013 

The authors have carried out MD simulations in lipid bilayer for three different models for the C-terminal 
region of the envelope glycoprotein E2. The use of different structures provides a priori a good 
benchmarking for the observations claimed in the article. However, the trajectories presented in this work 
seem not to have converged completely. The plots will show the trend more clearly if running average is 
used rather than the raw data. With the current advances in hardware and software, it seems to me that 
the authors have explored a limited conformational space of the system, which might be masking some 
other biological conformations for the region in the study. 

I am missing a more elaborated discussion in context with the available experimental data that they refer 
to in the conclusion. So far the correlation done in this study is too vague. The paper with experimental 
data the authors mentioned also suggests that this amphipathic helix could fold upon lipid binding. A MD 
simulation in water could help to elucidate this experimental observation. It would be convenient if the 
authors could specify the exact definition of membrane interface; for those people not working in the field, 
it could be a bit confusing. This concept could suggest the interface of the phospholipid bilayer 
(membrane core) or the lipid-water interface. 

The legend for the Supplementary Figure is missing the colour code used. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard, however I have significant reservations, as outlined 
above. 

Competing Interests: No competing interests were disclosed. 
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Martin Zacharias 

Physics Department, Technical University Munich, Garching bei Munchen, Germany 

Approved with reservations: 26 March 2013 
Referee Report: 26 March 2013 

The paper describes comparative MD simulations on the E2 TM region embedded in a lipid membrane. 
Different start conformations have been considered based on different molecular modelling methods + a 
structure from the pdb. I believe the simulations have been performed properly. Some interesting results 
on orientational preference of the helices have been reported. However, it should be emphasized that the 
simulations are too short to achieve convergence of the simulation results. This should be discussed in 
the Results & Discussion section or in the Conclusion section. Also, comparison to other simulation 
studies on peptide helices in membranes are missing (e.g. of the Thielemann group). At the end of the 
conclusion section the authors state "The results also revealed the orientation and positional preferences 
of the amphiphilic segment in relation to the water-bilayer interface that further clarify speculations from 
experimental studies." 

The authors mention ".. speculations from experimental studies.", however, no reference to such 
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speculations have been given and the details of the speculations are not discussed. The authors should 
give proper reference and should in detail explain what is meant by their statement. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard, however I have significant reservations, as outlined 
above. 



Competing Interests: No competing interests were disclosed. 
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